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Abstract. In this paper wc study planar polynomial differential systems of 
this form: 

— =X = A(X, Y), —=Y = B(X, Y), 
dt dt 

where A,B G Z[X,y] and degA < d, degB < d, \\A\\oc < H and \\B\\aa < 
"H. A lot of properties of planar polynomial differential systems are related 
to irreducible Darboux polynomials of the corresponding derivation: D = 
A{X,Y)dx + B{X,Y)dY ■ Darboux polynomials are usually computed with 
the method of undetermined coefficients. With this method we have to solve 
a polynomial system. We show that this approach can give rise to the compu- 
tation of an exponential number of reducible Darboux polynomials. Here wc 
show that the Lagutinskii-Pereira's algorithm computes irreducible Darboux 
polynomials with degree smaller than N, with a polynomial number, relatively 
to d, log('H) and N, binary operations. We also give a polynomial-time method 
to compute, if it exists, a rational first integral with bounded degree. 



Introduction 

In this paper wc study the foUowing planar polynomial differential system: 

^^X^A{X,Y), ^=Y = B{X,Y), 

where A,B e 1[X,Y] and deg^ < d, degS < d, \\A\\oc < H and ||B||oo < H. 
We associate to this polynomial differential system the polynomial derivation D = 
A{X, Y)dx+B{X, Y)dY. 

A polynomial / is said to be a Darboux polynomial, if D{f) = g./, where g is 
a polynomial. A lot of properties of a polynomial differential system are related 
to irreducible Darboux polynomials of the corresponding derivation D. Usually 
Darboux polynomials are computed with the method of undetermined coefficients. 
In other words, if we suppose that deg f < N then D{f) = g.f gives a polynomial 
system in the unknown coefficients of g and /. Then we can find / and g if we solve 
this system. We will see that this strategy can give rise to the computation of an 
exponential number of reducible Darboux polynomials. 

In this paper we show that we can compute all the irreducible Darboux polynomials 

,o(i)> 



of degree smaller than N with O \^[dN \og{'H)) j binary operations. Our strat- 
egy relies on the factorization of the ecstatic curve as suggested by J.V. Pereira in 
IPerOl) . 

This complexity result implies that, if we use the Prelle-Singer's strategy [PS83| . 



Date: September 16, 2010. 



1 



2 



G. CHEZE 



then we can compute an integrating factor in polynomial-time. With this integrat- 
ing factor we can then deduce an elementary solution of the given differential equa- 
tion. We also show that we can decide if D has a rational first integral of degree N. 
Furthermore, we can compute this rational first integral with o(^(^dNlog(l-L))'~'^^^^ 
binary operations. 



Related results. The computation of a first integral of a polynomial differential 
system is an old and classical problem. The situation is the following: we want to 
compute a function such that the curves J-{X,Y) = c, where c are constants, 
give orbits of the differential system. Thus we want to find a function T such that 
D{T) = 0. 

In 1878, G. Darboux |Dar78j gives a strategy to find first integrals. One of the tool 
developed by G. Darboux is now called Darboux polynomials. There exist a lot of 
different names in the literature for Darboux polynomials, for example we can find: 
special integrals, eigenpolynomials, algebraic invariant curves, particular algebraic 
solutions or special polynomials. 

Now recall briefiy the history and the use of Darboux polynomials. 
G. Darboux shows in |Dar78j that if we have enough Darboux polynomials then 
there exists a rational first integral for D, i.e. D{F) = and F G C(X, Y). More 
precisely, G. Darboux shows that if we have d{d -f l)/2 -I- 2 Darboux polynomials 
then the derivation has a rational first integral which can be expressed by means 
of these polynomials. This work is improved by H. Poincare jPoi91| . P. Painleve 
[Pai91j . and L. Autonne |Aut91j at the end of the XIX century. 

Since 1979, a lot of new results appear. In his book | Jou79| . J. -P. Jouanolou gives 
a polynomial derivation with no Darboux polynomials. Then there exist polynomial 
differential systems with no nontrivial rational first integral. 

In [PS83j M. Prelle and M. Singer give a structure theorem for polynomial differ- 
ential systems with an elementary first integral. Roughly speaking an elementary 
function is a function which can be written in terms of polynomials, logarithms, 
exponentials and algebraic extensions. With this structure theorem the authors 
show how we can compute an elementary first integral if we know all the Darboux 
polynomials of D. 

In |Man93j and jMM97j Y.-K Man and M. MacCallum explain how we can use 
Prelle-Singer's method in practice. The bottleneck of this method is the computa- 
tion of Darboux polynomials. Indeed, the method of undetermined coefficients is 
used. With this method we have to solve a polynomial system, and Y.-K Man in 
[Man93| explains how to solve this system with a Groebner basis. The resolution 
of this polynomial system is difficult because no particular structure of this system 
is known. Furthermore, this polynomial system can give an exponential number of 
reducible Darboux polynomials. 

In [Sin92) M. Singer shows a structure theorem for polynomial differential sys- 
tems with a Liouvillian first integral. Roughly speaking a Liouvillian function 
is a function which can be obtained "by quadratures" of elementary functions. 
C. Christopher in jChr99| improves this structure theorem, and then he suggests 
an algorithm to find Liouvillian first integrals. 
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In a series of papers jDDdMS02| IDDdM02a| IDDdM02bj the authors describe an 
algorithm and its implementation for the computation of Liouvillian first integrals. 
The strategy is in the same spirit as the one proposed by C. Christopher. 

Singer's theorem, Christopher's theorem, and Duarte-Duarte-da Mota's algo- 
rithm are based on Darboux polynomials. 

In |CLP07| the authors define an algebraic, a geometric, an integral, an infini- 
tesimal and an holonomic notion of multiplicity for an irreducible Darboux poly- 
nomial. They show that these notions are related. Furthermore, in order to define 
the algebraic multiplicity they introduce the ecstatic curve and then use some of 
its properties. We recall the definition and some properties of the ecstatic curve in 
Section [3l It seems that this curve and some of its properties was already known by 
M.N. Lagutinskii (1871-1915), see jDLS98| . and was rediscovered by J.V. Pereira, 
see [PerOlj . Theorem 2 in |Per01| is the main tool of our algorithms. 

In [CMS091 ICMSOGj , the authors give a Las Vegas strategy to decide if there 
exist Darboux polynomials for a given derivation. 

Recently, in jFGlO] the authors propose an algorithm to compute a rational 
first integral without computing Darboux polynomials. Unfortunately, there are no 
complexity study for this algorithm. 

Darboux polynomials are also used in the qualitative study of polynomial system. 
For example, the inverse integrating factor is a special Darboux polynomial, and, 
algebraic limit cycles are factors of the inverse integrating factor, see e.g. jGLV96j . 

We also mention that Darboux polynomials are used in Physics: In |Hew91j the 
author uses Darboux polynomials in order to find exact solutions to the Einstein 
field. In (LZOO, Val05, TV08. Darboux polynomials are also used to study the Rik- 
itake system, which is a simple model for the earth's magneto-hydrodynamic. 

In |CG03| IGinOTj the reader can find some open questions and some relations 
between the computation of Darboux polynomials and Hilbert's 16th problem. 

For other results, references or applications of Darboux polynomials the reader 
can consult for example jGorOll IDLA06| . 

Main results. In this paper we show: 

Theorem 1. Let D = A(X,Y)dx + B{X,Y)dY be a polynomial derivation such 



that A{X,Y),B{X,Y) e Z[X,Y], degA < d, degB < d, \\A\\^ < H, \\B\\^ < % 



and A, B are coprime. 

(1) We can decide if there exists a finite number of irreducible Darboux poly- 
nomials with degree smaller than N in a deterministic way with 



(2) If there exists a finite number of irreducible Darboux polynomials with degree 
smaller than N then we can compute all of them in a deterministic way with 



To author's knowledge this is the first polynomial-time result on the computa- 
tion of Darboux polynomials. 





4 



G. CHEZE 



In this paper we suppose that A and B have integer coefficients and are coprime. 
This hypothesis is not restrictive. Indeed, if A and B have rational coefficients then 
we can always reduce our study to the case of integral coefficients by clearing de- 
nominators. Furthermore if A and B have a nontrivial greatest common divisor 
R then wc set ds = Rdt. This gives a derivation D2 = A/Rdx + B/Rdy- If D2 
has a first integral then D has a first integral, thus the hypothesis "A and B are 
coprime" is not restrictive. 

Furthermore in this paper we describe algorithms in the bivariate case in order 
to emphasize their role in the Prelle-Singer's algorithm. Nevertheless, Proposition 
[18] which is our main tool is also true in the multivariate case, i.e. when we con- 
sider a derivation D = Aidxi + • • • + Andx„ where Ai £ Z[Xi, . . . , Xn]- Then the 
Lagutinskii-Pereira's algorithm, see Section |4l is also correct in the multivariate 
case. 

If the polynomial derivation D has an infinite number of Darboux polynomials 
then by Darboux's theorem, D has a rational first integral. Thus in this situation 
the problem is the computation of a rational first integral. We prove the following 
result: 

Theorem 2. Let D = A{X,Y)dx + B{X,Y)dY be a polynomial derivation such 
that A{X,Y),B{X,Y) e Z[X,Y], degA < d, degB < d, \\A\\oo < U, \\B\\oo < H 
and A, B are coprime. 

(1) We can decide if there exists a rational first integral J- of degree smaller 
than N with ©^(dA^ log('H))'^''^'^ binary operations. 

(2) // there exists a rational first integral with degree smaller than N then we 
can compute it in a deterministic way with o((d7Vlog(H))'^^^^) binary 
operations. 

To author's knowledge this is the first polynomial-time result on the computa- 
tion of rational first integrals. 



Structure of this paper. In Section [T] we recall some classical results about 
Darboux polynomials, the spectrum of a rational function and the complexity of 
bivariate factorization. In Section [2] we show that the method of undetermined 
coefficients can give an exponential number of reducible Darboux polynomials. In 
Section[3]we give the definition and some properties of the ecstatic curve. In Section 
m we prove Theorem [T] and in Section O we prove Theorem [21 At last in Section [5] 
we ask two questions about complexity and polynomial differential equations. 

Notations. Let f{X, Y) = - f,,jX'Y^ e Z[X, Y] be a polynomial. 
||/||co = i^f^Xijlfijl is the height of the polynomial /. 

degf is the total degree of the polynomial /. The degree of a reduced rational 

function p/q is the maximum of degp and degq. 

The bit-size of a bivariate polynomial / is (deg /)^ log(||/||oo). 

dx (resp. dy) denotes the derivative relatively to the variable X (resp. Y). 

We set: div(^, B) = dxA + dyB. 
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1. Classical results 

1.1. Darboux polynomials. 

Definition 3. Let D — A{X,Y)dx + B{X,Y)dY be the derivation associated to 
the planar polynomial differential system 

X = A{X,Y), Y^B{XX). 

where A,B 7L\X, Y\. 

A polynomial / e C[X, Y\ is said to be a Darboux polynomial of D if there exists 
a polynomial g E C[X, Y] such that D{f) = g.f. 
The polynomial g is called the cofactor of /. 

If is a function such that D{F) = then we said that is a first integral of D. 

Proposition 4. Let f he a polynomial and let / = /1/2 be a factorization of f 
where fi and /2 are coprime. Then f is a Darboux polynomial with cofactor g if and 
only if fi and /2 are Darboux polynomials with cof actors gi and g2 ■ Furthermore 
9 = 91 +92- 

Proof. Sec for example Lemma 8.3 page 216 in [DLA06j . □ 
Now we recall Darboux's Theorem. 

Theorem 5 (Darboux's Theorem). Let A,B e Z[X, Y] and let D = Adx + Bdy- 

If fiT ■ ■ T fm G C[X, y] are relatively prime irreducible Darboux polynomials for 
i = 1, . . . , m, then either m < d(d + l)/2 + 2 where d = max(deg A, deg B) or there 
exist integers Hi not all zero such that D{w) — 0, where w — YYiLi fj^' ■ In the latter 
case, if f is any irreducible Darboux polynomial, then either there exists A, jj, in <C, 
not both zero such that f divides ^Yiiei fi^ ~ f-Wj^j fj where / = {i | rii > 0} 
and J ~ {j \ Uj < 0}, or G divides gcd(^, B). 

Proof. See |Dar78| or |Sin92| . □ 

Definition 6. A function R is an integrating factor if 

D{R) = -div{A,B)R. 

Remark 7. We remark that an integrating factor satisfies also one of these equivalent 
conditions: 

diY{RA,RB) = 0, dx{RA) = -dyiRB). 

If we know an integrating factor then we can deduce a first integral. Indeed, if 
R is an integrating factor then 

(1.1) ^ J [rA + Oy j RBdX^dY 

is a first integral. 

1.2. Elementary solutions. The following theorem is due to M. Prelle and M. Singer, 
see [PS83] . 

Theorem 8 (Prelle-Singer). If a polynomial differential system has an elementary 
first integral then there exists an integrating factor which is a Kth root of a rational 
function. 
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This gives the foUowing method to compute an elementary first integral. This 
method is called the Prelle-Singer's method. Here we follow the description given 
in |Man93j . 

Prelle-Singer's method 

Input: D = A{X^ Y)dx + B{X, Y)dY a polynomial derivation, and N an integer. 
Output: An elementary first integral, if it exists, constructed with Darboux poly- 
nomials of degree smaller than N . 

(1) Set n=l. 

(2) Find all monic irreducible polynomials fi such that deg fi < n and fi divides 
D[h). 

(3) Let D{fi) = Qifi. Decide if there are constants n^, not all zero, such that 

m 

^ rug^ = 0. 
1=1 

If such rii exist then n"=i * integral. 

If no such rii exist then go to the next step. 

(4) Decide if there are constants n^, such that 

m 

'^rngi = -div(yl, B). 

i=l 

If such Ui exist then J^™ is an integrating factor for the given differ- 
ential equation. A first integral is given by the formula (jl.ip . 
If no such rii exists then go to the next step. 

(5) Increase the value n by 1. If n is greater than then return failure other- 
wise repeat the whole procedure. 

If we want to compute an integrating factor, it is well known, see |Man93| . that 
step m is the most difficult step of the Prelle-Singer's method. Indeed, when we 
have all the irreducible Darboux polynomials we just have to solve linear systems 
(see Step O and Step [3]) to deduce and integrating factor. 

Remark 9. In Prelle-Singer's method the user must give a bound A'^. Nowadays we 
cannot remove this input. Indeed, we do not know a bound on the maximal degree 
of irreducible Darboux polynomials of a given derivation. This is an open question 
and appears in |Poi91[ IPS83| . The following example shows that we cannot get a 
bound in term of the degree of D only. The bound must also take into account the 
coefficients of D. 

Example: The derivation D = {n + l)Xdx + nYdy has A" — y+i as Darboux 
polynomial. 

1.3. Spectrum of a rational function. In this subsection we recall the definition 
and a property of the spectrum of a rational function. This notion is used in Section 

m 

Definition 10. A rational function f{X,Y) e Q(X,Y) is said to be composite if 
it can be written f = uo h where h{X,Y) e Q{X,Y) and u € Q(T) such that 
deg(u) > 2. Otherwise / is said to be non-composite. 
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Definition 11. Let f = p/q ^ Q{X,Y) be a reduced rational function of degree 

d. The set 

a{p, g) = {(A : Ai) e Pc | \p + nq is reducible in C[X, Y], 

or deg(Ap + fiq) < d} 

is the spectrum of / = p/q. We recall that a polynomial reducible in C[X, F] is 
said to be absolutely reducible. 

The spectrum cr(p, q) is finite if and only ii p/q is non-composite and if and only 
if the pencil of algebraic curves Xp + fiq = 0, has an irreducible general element 
(see for instance |Jou79| Chapitre 2, Theoreme 3.4.6] or jBod08[ Theorem 2.2] for 
detailed proofs). 

To author's knowledge, the first effective result on the spectrum has been given 
by Poincare |Poi91j . In this paper, Poincarc gives a relation between the number 
of saddles of a polynomial vector field and the spectrum. He also shows that 
WiP, q)\ ^ (2rf — 1)^ + 2d + 2. This bound was improved only recently by Ruppcrt 
|Rup86| who proved the following result: 

Proposition 12. If p/q £ Q{X,Y) is a reduced non-composite rational function 
of degree d then \a{p,q)\ < d^ — 1. 

This result was obtained as a byproduct of a very interesting technique developed 
to decide the reducibility of an algebraic plane curve. 

Several papers improve this resuh, see e.g. |Lor93[ IVis93l IAHS03[ [BodOSl IBC08| . 

Remark 13. If p/q is a reduced non-composite rational first integral of D then p+Xq, 
where A € C, are Darboux polynomials. Then D has infinitely many irreducible 
Darboux polynomials because the spectrum a{p, q) is finite. 

1.4. Complexity results. In this paper we consider the dense representation of 
polynomials in the usual monomial basis. We recall that we can factorize in a de- 
terministic way an integer univariate polynomial / with 0^((ilog('H))'^''^'^ binary 

operations, where d is the degree of / and H its height, see e.g. |Sch84j . An algo- 
rithm with this kind of complexity is called a polynomial-time algorithm because 
its complexity is bounded by a polynomial in the size of the input. We recall that 
the bit-size of an univariate polynomial / of degree d and height "H is d\og{T-L). 
The first polynomial-time algorithm for the factorization of univariate polynomials 
is due to Lenstra, Lenstra and Lovasz, see |LLL82| . However this algorithm, called 
LLL, is probabilistic because it uses Berlekamp's algorithm which is probabilistic, 
see |Ber70j . Nevertheless, there exist deterministic polynomial-time algorithms, see 

e. g. |Sch841 |KLL88| . In jSch841 lKLL88j the strategy is numerical. Instead of com- 
puting a modular factorization as in jLLL82| . the authors propose to compute a 
complex root with a sufficiently high precision. 

In |Kal85bj the author shows that we can reduce bivariate factorization to uni- 
variate factorization. Then we get a deterministic polynomial-time algorithm for 
the factorization of integer bivariate polynomials. Another polynomial-time algo- 
rithm is proposed in jLen84| . the authors extend to the multivariate case the LLL 
algorithm. 

Few time later, several papers see |Kal85a| IDT89| ITra85| . show that we can also 
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compute the absolute factorizatiou (i.e. the factorization in Q[X, F], where Q is 
the algebraic closure of Q) of integer bivariate polynomials with a deterministic 
polynomial-time algorithm. Thus there exist deterministic algorithms which per- 
form the absolute factorization of a bivariate polynomial f{X,Y) € Q[X, y] of 

degree d and height H with o(^{dlog{'H))'^'^^^^ binary operations. 
For the history of factorization's algorithm the reader can consult |Kal90[ IKal92j . 
For more recent results about complexity and bivariate polynomials factorization 
see |Gao03l IBLS+Ol ILecOGl IB^KSOQ. ,CL07l IAKS071 . 

We also mention here that we can solve a linear system, and compute a deter- 
minant, with coefficients in Z or in Z[X, y], in polynomial-time, see e.g. |YapOO[ 
IBCS97[ lHS75j . This will be useful in our complexity analysis. 

2. The method of undetermined coefficients 

Usually, Darboux polynomials are computed with the method of undetermined 
coefficients. In other words, we write: D{f) = g.f with degg < d, where d is 
the degree of D, and deg/ < N. This gives a polynomial system in the unknown 
coefficients of g and /. This system has 0{N'^ + d^) unknowns. In this section, 
we show that the method of undetermined coefficients can give an exponential 
numbers of reducible Darboux polynomials. We recall that only irreducible Darboux 
polynomials are useful because the product of Darboux polynomials is a Darboux 
polynomial, see Proposition |4l 

Lemma 14. // we consider the following derivation: 

d-l 

D = {dYJ')dx - {dxT)dY, with T{X, F) = F J| (X + i) + X, 

1=1 

then there are at least 2'^^^ + 1 Darboux polynomials with degree smaller than d. 

Proof. In this situation X + i, i = 1, . . . , d— 1, are irreducible Darboux polynomials. 
Thus 

+ i), where / is a subset of {1, . . . ,d~ 1} 

I 

is a Darboux polynomial. We have 2'^~^ such Darboux polynomials and these 
polynomials have a degree smaller than d — 1. Furthermore J-' is an irreducible 
Darboux polynomial and deg\F = d. This gives the desired result. □ 

Remark 15. In Lemma 1141 we give a derivation with a rational first integral. It 
would be interesting to have the same kind of result with a derivation with no 
rational first integral. 

In the situation of Lemma [T?l J- is an irreducible Darboux polynomial of degree 
d. Thus if we want to find all the irreducible Darboux polynomials, then we have to 
consider polynomials / with degree smaller or equal to d in the polynomial system 
D{f) = g.f. Now, Lemma |14I implies that we have at least 2"^^^ + 1 solutions for 
the system D{f) = g.f with deg f < d. 

Then we can conclude: in the worst case the method of undetermined coefficients 
gives an exponential number, in d, of reducible Darboux polynomials. 
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The problem comes from reducible Darboux polynomials. The recombination of 
irreducible Darboux polynomials gives an exponential number of Darboux polyno- 
mials. We can "avoid" this problem. Indeed, we can add to the system D{f) = g.f, 
forms ipi,iq{f), i = 1, . . . , T, such that fi,N{f) = 0, for all i = 1, . . . , T if and only 
if / is an absolutely reducible polynomials with degree N. Such forms exist and 
they are called Noether's forms, see |Rup86[ ISchOOj . Unfortunately, to author's 
knowledge, the number T of these forms is exponential in N. More precisely, the 
number of Noether's forms is equal to ( jv^-i }■ With Stirling's formula we can 
show, when N tends to infinity: 

/2N^ -3N + 1\ 2^^'-3^+i 

Thus in this case we have an exponential number of equations, and then we get a 
method with an exponential complexity in N. 

We remark that there exists a strategy to compute the leading term of Darboux 
polynomials, see e.g. |Chr94j . This strategy gives a net gain in practical examples, 
see |MM97| . However, this strategy do not detect the leading term of an irreducible 
Darboux polynomial. Then, even if we add this strategy to the method of undeter- 
mined coefficient, we still get an exponential number of Darboux polynomials. 

The exponential complexity is related to the recombination of irreducible Dar- 
boux polynomials. Exponential complexity due to recombinations appears also 
when we study factorization algorithms. However, we can factorize polynomials in 
polynomial-time as mentioned in Section [1.41 In the following we show that we can 
reduce the computation of irreducible Darboux polynomials to the factorization of 
a bivariate polynomial. Then, we will deduce a polynomial-time algorithm for the 
computation of irreducible Darboux polynomials. 



3. The ecstatic curve 



Definition 16. Let D he a. polynomial derivation, the iVth ecstatic curve of D, 
^b,n{D), is given by the polynomial 



det 



Vl 

D{vi) 



V2 
D{V2) 



D{vi) 



\ 



where S = {wi , W2 , ■ . ■ , } is a basis of C [X, F] <Ar , the C- vector space of polynomials 
mC[X,Y] of degree at most N, I = {N + 1){N + 2)/2, and D''{v,) = Z)(D'=-i (w;)) ■ 
When B is the monomial basis, we denote by £n{D) the ecstatic curve. 

Now, we remark that the ecstatic curve is independent of the chosen basis of 
C[X,Y]<:pf up to a multiplicative constant. 

Proposition 17. Let B be a basis of C[X,Y]<n . We have 
where c G C* . 
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Proof. We are going to show that Eb,n{D) is the determinant of an endomorphism. 
Consider 

V:<C[X,Y)[U,V]<N <C{X,Y)[U,V]<N 

U^yi I y ^ £)cantor(fe,/)^j^iyj^jjfc^/ 

Q<k+l<N 

where U and V are new independent variables, C{X,Y)\U,V]<n is the C(X, F) 
vector space of polynomials with coefficients in C(X, F), and with a degree relatively 
to U and V at most N . Furthermore cantor is the following map: 

cantor : — > N 

[kj) I — > . 

This application maps (0, 0) to 0, (0, 1) to 1, (1, 0) to 2, (0, 2) to 3. . . That is to 
say, the monomials U'^V'' are ordered with a graded lexicographic order. 

2? is defined on a basis thus I? is a well defined endomorphism. £n{D) is the 
determinant of this endomorphism written in the monomial basis. 

Now we remark that: 

0<i+j<N 0<k+l<N 0<i+j<N 

where pij € C. 

Each element of the basis B can be written j Pi.j^^^'' with j E C. We can also 
consider a basis B' of the C{X,Y) vector-space C{X,Y)[U,V]<n, such that each 
element is written J^i j Pi,jU^V^ with pi_j £ C. Thus £b,n{D) is the determinant of 
2? written with B' in the domain, and with the monomial basis in the target space. 
As Pi J e C, this gives the desired result. □ 

The following proposition is due to J.-V. Pereira, see [PerOll Proposition 1]. 
It is the key point of our algorithm: it shows that the computation of Darboux 
polynomials can be reduced to the factorization of £n{D)- We give a proof in order 
to ease the readability of the paper. 

Proposition 18. Every Darboux polynomial, relatively to D, of degree smaller 
than N is a factor of £n{D). 

Proof. Let F S C[X,Y]<i\[ be a Darboux polynomial. By Proposition [T71 we can 
choose a basis B where vi = F. 
Furthermore, we have: 

D{F) = giF, 

D\F) = D{cnF) = {gl+D{cn))F = g2F, 



D'-\F) = .g;_iF, 
where gi, g2, ■ . ■ , gi-i arc polynomials. 

Thus F is a factor of £b.n{D) and this concludes the proof. □ 

Remark 19. The converse is false. Indeed, consider the derivation D = —2X^dx + 
(1 - AXY)dY. Then £i{D) = YX* but Y is not a Darboux polynomial. 
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We know by Darboux's Theorem, see Theorem \E\ that if a derivation has a ra- 
tional first integral then there are infinitely many irreducible Darboux polynomials. 
Thus if D has a rational first integral then by Proposition [T51 £n{D) has infinitely 
many irreducible factors. This gives £n{D) ~ 0, and also £m{D) = for M bigger 
than N. The following proposition says that the converse is also true. This propo- 
sition will be useful in Section [5] when we will study the computation of rational 
first integrals. 

Proposition 20. We have £iq{D) = and £n-i{D) ^ if and only if D admits 
a rational first integral of exact degree N . 

Proof See |Per01[ Theorem 1]. □ 

In our complexity study we will need to know the bit-size of £n{D). Wc recall 
that the bit-size of £Ar(i:>) is: {dcg{£N{D)) log (||£Ar((i)||oo). Thus in the following 
we are going to compute the degree and the height of £n{D). 

Proposition 21. Let D = A{X,Y)dx + B{X,Y)dY he a polynomial derivation, 
where degA,dcgJ5 < d. The degree of £n{D) is at most N.l + {d — 1).{1 — l).l/2, 
where I = {N + 1){N + 2)/2. 

Proof. By definition of a determinant we have: 

i~i 

degf^p) <^deg 

fc=0 

where D^(vi) = Vi. 

A straightforward computation gives AcgD^{vi) < k{d- l) + N. Then 
degfjvP) < ^(fc(d-l) + 7V), 

fe=0 



< N.l + {d-l)J2k, 



k=0 

< N.l + {d-l).{l-l).l/2. 
This gives the desired result. □ 

Corollary 22. Under the hypothesis of Provosition [2l[ we have deg£N{D) belongs 
to 0{dN^). 

Proposition 23. The height |15jv(£')||oo satisfies 

OS i(i-i)/2 



\\£n{D)\\oo < (2lH{l{d-l)+Nf) 



Proof. First wc recall that if /i and /2 arc two polynomials with total degree smaller 
than d then we have 

|!/l./2||oo<(d+l)'|l/l|loo||/2||oo. 

This gives for / e C[X,Y]<n, 

\\D{f)\\^<2nN'\\f\U 
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By induction, using {N + k{d - 1) + 1)^ < (iV + k{d - l))^, we get 

i-l 

\\D\f)\\oo < {k{d - 1) + Nf) ll/lloo, where i > 1. 



fc=0 

By definition of a determinant, we have: 

i-i 

\\£NiD)\\^ < j2 nii^'(M^+i))ii-' 

(tGSi 1=0 

where {vi, . . . ,vi} is the monomial basis of C[X, Y]<:n. We get then: 

(-1 i-l 

\\Sn{D)\\^ < l\l[TW(^l[{k{d-l) + Ny 

i=l k=0 

i-i 

3i 



< l\Y[2'Wi^{l{d-l)+N 

i=l 

/ q\'('-l)/2 

< [2lH{l{d - I) + N) j 



□ 



Corollary 24. Under the hypothesis of Proposition [23[ the bit-size of £n(D) be- 
longs to o(^{dN\og{n))°''^^y 

Proof. By Corollary [22l and Proposition [23l the bit-size of £n{D) is bounded by 

{dN^f.l.{l - 1). log (2m{l{d - 1) + iV)^) . 
This gives the desired result. □ 

4. D HAS A FINITE NUMBER OF IRREDUCIBLE DARBOUX POLYNOMIALS 

If D has a finite number of irreducible Darboux polynomials then the ecstatic 
curve £m{D) is non-zero by Proposition 1201 It follows that we can compute ir- 
reducible Darboux polynomials with a bivariate factorization algorithm thanks to 
Proposition [iHl It seems that Proposition [TSl has been proved by M.N. Lagutinskii 
and rediscovered by J.V. Pcreira, see [PerOlj Theorem 2]. Then we call "Lagutinskii- 
Pereira's algorithm" the following algorithm. 

Lagutinskii-Pereira's algorithm 

Input: A polynomial derivation D = A{X,Y)dx + B{X,Y)dY , and N an integer. 
Output: The finite set S of all the absolute irreducible Darboux polynomials with 
degree smaller than A'^ or "There exists an infinite number of irreducible Darboux 
polynomials" . 



(1) S = {). 

(2) Compute £n{D). 

(3) If £n{D) — then Return "There exists an infinite number of irreducible 
Darboux polynomials" else go to step 21 end If. 

(4) Compute the set /i, . . . , of all absolutely irreducible factors of £n{D) 
with degree smaller than N . 
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(5) For z := 1, . . . , m do: If gcd(/i, D{fi)) = fi then add fi to S, end If, end 
For. 

(6) Return 5". 

Proposition 25. The Lagutinskii-Pereira's algorithm is correct. 

Proof. This is a straightforward consequence of Proposition [181 O 

Now we can prove Theorem [T] 

Proof. The Lagutinskii-Pereira's algorithm is correct and works with the claimed 
complexity. Indeed, we can compute £n{D) in polynomial-time because it is a 
determinant. Furthermore, by Corollary [24] we know that the bit-size of £n{D) 

belongs to 0^(diVlog('H))*^'^''y Then we can compute the absolutely irreducible 

factors of £n{D) with 0^(rfiVlog('H))'^^^''^ binary operations, see Section \TA\ As 

gcd computations can also be performed in polynomial-time we obtain the desired 
result. □ 

We also deduce: 

Corollary 26. Under conditions of Theorem]^ if there exists an integrating factor 
R such that R = J^^ , where fi are Darboux polynomials with degree smaller 

than N , then we can compute R with 0^(rfiVlog('H))'^''^''^ binary operations. 

Proof. We compute Darboux polynomials with the Lagutinskii-Pereira's algorithm 
and then we solve a linear system as in Step |4| of the Prclle-Singer's method, see 
Section [H □ 

This corollary implies that we can compute an integrating factor corresponding to 
an elementary first integral with the Prelle-Singer's method with O ^ {dN log('H)) '^^^^ 
binary operations. 

5. D HAS A RATIONAL FIRST INTEGRAL 

We have seen in Section [3| that is easy to test if a derivation has a rational first 
integral with degree smaller than N . Indeed, we just have to compute the ecstatic 
curve £n{D) and check if it is the zero polynomial. In this section, we show how 
we can compute this rational first integral. 

5.1. Computation of a Darboux polynomials with degree N. In this section 
we suppose that D has a rational first integral p/q with degree N. Then p and q 
are Darboux polynomials with degree N. By Proposition [201 £n{D) = 0. Thus we 
cannot compute p and q as factors of £n{D)- The strategy is then the following: 
compute one Darboux polynomial with degree N , compute its cofactor, and then 
deduce p and q. 

Now we explain how we can compute a Darboux polynomial with degree N . 

Definition 27. Let D be a polynomial derivation, £^^{0), is the polynomial 
^b°,n{D) where is the monomial basis of the C- vector space of polynomials in 
C[X, y] of degree at most with constant term equal to zero. 

We have the following property. 
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Proposition 28. Letp/q be a non-composite rational first integral of D, such that 
degip/q) = N, and p{0, 0)g(0, 0) 7^ 0. 

(1) We have £n,o{D) ^ in Q[X, Y]. 

(2) If we set (Ao,/io) = ( — 9(0j 0),p(0, 0)) then Xap+fioq is a factor of £n,o{D). 

Proof. This proof follows very closely the proof of Theorem 5.3 and Proposition 5.2 
in |CLP07j . 

First, we prove that £n,o{D) ^ 0. 

We suppose the converse: £n,o{D) = and we show that it is absurd. 

If £n.o{D) ~ then the columns of the matrix are linearly dependent. Hence there 

are rational functions Ci{X,Y) € Q{X,Y) such that 

k 

(5.1) Af, c^DHv^) = 0, J = 0, . . . , fc - 1 

1=1 

with k = {N + 1){N + 2)/2 — 1. Now, take k to be the smallest value such that for 
i = 1, . . . ,fc there exists rational functions Ci, not all zero, and w,; € C[X, y]<jv,Oi 
linearly independent over Q such that equalities (|5.ip holds. It is clear that k > 1. 
We have: 

k 

Z?(AA,) - A/;+i = ^ Z?(c,)i?n"«) = 0, J = 0, . . . , fc - 2, 

i=l 

and so from the minimality of fc, we see that the terms D(ci) must all vanish. 
Hence, each of the Ci are cither rational first integrals or constants. 
Now, we consider the polynomial G = Xp + fiq where (A : fi) does not belongs to 
(t(p, g) U {(Ao : A*o)}- This choice is possible because by Proposition [121 <^{P:Q) 
is finite. Then G is absolutely irreducible, deg G = N and G(0, 0) ^ 0. We also 
remark that Ci are constants for all (x, y) such that G{x,y) = 0, because p/q is a 
first integral. We denote by Ci{X,fj.) these constants. This gives 
fc 

Nq''^ := ^ c,;(A, fi)vi 7^ 0, and G = Xp + ^iq divides Mq'^ . 

1=1 

Indeed, Vi are linearly independent, and furthermore if G{x, y) = then A/j^''^ (x, y) = 
No{x, y) = 0. Thus, there exists c € Q such that c.G = A/jf'^. As by construction, 
<N,Q, we deduce that G belongs to C[X, F]<jv,o and this contradicts 
G(0,0) 7^ 0. Thus we obtain EnA^) ^ 0. 

Second, we prove that F = X^p + p^^q divides £n o{D). 

is a Darboux polynomial and belongs to C[X, y]<jv,o- Since £n,o{D) is indepen- 
dent of the chosen basis (with the same arguments used in Proposition [17)) . we can 
choose a basis where vi = F. Then as in Proposition [T8| we deduce that F is a 
factor of £]y q{D) and this concludes the proof. □ 

5.2. Computation of a rational first integral. Thanks to Proposition [28l we 
can describe an algorithm which computes a rational first integral. 

We denote by {xk,yk), fc = 1, . . . , N^, the points in S* x 5 where 
S = {0, . . . , — 1} and by Dk the following derivation: 

Dk = A{X + Xk,Y + yk)dx + B{X + Xk, Y + yk)dY. 
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Let g e Q[X, y] with degree smaller than d, Ck.g is the following linear map: 

f ^ Dk{f)-g.f 

Lemma 29. Suppose that Dk has a rational first integral p/ q of degree N . We 
denote by g the cof actor of p and q. 

Then dimQker£fc^g = 2 and if we denote by {p^q} a basis ofkeiCk^g thenp/q is a 
rational first integral of ■ 

Proof. This proof follows very closely the first part of the proof of Theorem 6 in 
|MO04! . 

If / G ker£fc.g then by Darboux's theorem, see Theorem [51 there exist a and (3 in 
Q such that f = ap + f3q or f has degree less than N and divides ap + Pq. The 
last case implies that (a,/?) S a{p,q). We deduce that: 

(5.2) ker Ck,g = Span{p,q) U Span{fi) U • • • U Span{fm), 

where fi is a factor of ap + (3q, {a, 13) G (y{p, q), and the cofactor of fi is equal to g. 
Remark that this union is finite because cr(p, q) is finite. 

As Q is infinite, equality (|5.2p implies "kcv Ck,g = Span(p,q). □ 

Algorithm Rat-First-Int 

Input: A polynomial derivation D = A{X, y)dx + B{X, F)9y , and N an integer. 
Output: A rational first integral with degree smaller than TV or "There exists no 
rational first integral with degree smaller than TV" . 

(1) Compute En{D). 

(2) If En{D) 7^ then Return "There exists no rational first integral with 
degree smaller than iV" , else go to step [3l end If. 

(3) Compute the smallest integer n such that £n{D) = and £„-i(D) ^ 0. 

(4) Set F:=0; k:=l; 

(5) While F = do 

(a) Compute £'„,o(£'fc)- 

(b) Compute all the irreducible factors /i, ...,/„ of £nfi{Dk) with degree 
equal to n. 

(c) For all i:=l, . . . , m do: If gcd (/i, Dk{fi)) = fi then set F := fi and 
go to step [51 end If; end For. 

(d) k:=k + l; 
end While. 

(6) Compute the cofactor g := Dk-i{F)/F. 

(7) Compute a basis {p,q} of ker£fc_i^g. 

(8) Return p/q{X - Xk-i,Y ~ yk-i). 

Proposition 30. The algorithm Rat-First-Int terminates and uses the While loop 
at most times. Furthermore the algorithm Rat-First-Int is correct. 

Proof. The algorithm terminates. We just have to show that the While loop termi- 
nates. 

In Step[5l thanks to Proposition[20l D has a reduced rational first integral p/g of de- 
gree n < N. Then Dk has a rational first integral /(/^ (AT, Y) = p/q{X-\-Xk, Y-\-yk) 
of degree n < N. 
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Furthermore p/q is non-composite (thus pk/qk is also non-composite). Indeed, if 
p/q is composite then p/q ^ u o h with deg h < n. This gives 

= D{p/q) = D{u oh)= u'{h)D{h), 

and u'ih) ^ because degu > 2. Therefore D{h) = and ft, is a rational first 
integral with degree smaller than n. Thus by Proposition [20l £dcgh{D) = 0, this 
contradicts the minimality of n. We deduce then: p/q is non-composite. 

Now, remark that the algorithm terminates if the following two conditions are 
satisfied: 

(5.3) p{xk,yk) ^ or q{xk,yk) ^ 0, 
and 

(5.4) ( - q{xk,yk) ■■ p{xk, Vk)) ^ cr{Pk,qk) cr{p, q)- 

Indeed, in this situation we can apply Proposition [25] and we deduce that there 
exists a polynomial F = -q{xk,yk)Pk{X,Y) + p{xk,yk)qk{X,Y) absolutely irre- 
ducible such that gcd{F, Dk{F)) = F. 

Now we show that there exists a point {xk,yk) in {0, . . . — 1}^ such that 
()5.3p and (|5.4p are satisfied. 

By Bezout's Theorem we just have to avoid N'^ points to satisfy ()5.3|) . 
Now, we consider the polynomial 

r{x,Y)^ J] {\p{x,Y) + ^Iq{x,Y)). 

(A:^i)e(T(p,g) 

We remark that if V{xk,yk) 7^ then (|5.4p is satisfied. Furthermore, by Proposi- 
tion [H degT' < N{N'^ - 1). Zippel-Schartz's lemma, see jvzGG03[ Lemma 6.44], 
implies that V has at most - roots in {0, . . . , A^^ - 1}^. Then the algorithm 
terminates and uses the While loop at most — + N"^ + I times. 

The algorithm is correct. If fi satisfies gcd (/j, Dk{fij) = /; then fi is a Darboux 
polynomial. Then by Darboux's theorem, see Proposition [SJ we deduce that fi = 
apk + Pqk because dcg pk/qk = n = dcg/i. 

As Pk/qk is a first integral, the cofactors of pk and qk are equal. Then, we deduce 
that the cofactors of pk, qk and fi are equal. This cofactor is the polynomial g. 
Then thanks to Lemma [29] we conclude that the algorithm is correct. □ 

Now we can prove Theorem [2l 

Proof. Thanks to Proposition [30] we just have to prove that the algorithm Rat- 
First-lnt works with O (^(dN log{'H))'~"^^^^ binary operations. 

We have already mention that we can compute determinants and solve linear sys- 
tems in polynomial-time. Thus we can perform Step [1] Step [2] and Step [3] with 

C'((rfiVlog(-H))*^^^^) binary operations. 

Now, we study the While loop. 

We recall here that we can shift the variable of a polynomial in polynomial-time see 
e.g. jBP94[ Problem 2.6]. We can also perform linear algebra, compute gcd and di- 
vide polynomials in polynomial-time, see e.g. |vzGG03] . Furthermore, as in Corol- 
lary [24] we can show that the bit-size of £jv,o(^) belongs to O ( (^dN \og{'H))'^^^^j . 
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Then we can factorize £n,o (D) with o(^{dNlog{n))°^^^^ binary operations. As 
we use the While loop at most times we obtain the desired result. □ 

6. Open questions 

6.1. Liouvillian first integrals. In this paper we have shown how to compute effi- 
ciently Darboux polynomials. Thus this improves the complexity of Prelle-Singer's 
method. Now the question is; Can we compute efficiently an integrating factor 
corresponding to a Liouvillian first integral? 

In |DDdMS021 IDDdM02a[ IDDdM02b| the authors give an algorithm to compute 
such integrating factors. The key point is the computation of exponential factors. 
The definition of an exponential factor is the following: 

Definition 31. Given f,g £ C[X,Y], we say that e = cxp{g/f) is an exponential 
factor of the derivation D if D{e)/e is a polynomial of degree at most d ~ 1. 

In [CLP07) the authors define the integrable multiplicity and the algebraic mul- 
tiplicity. 

Definition 32. We say that a Darboux polynomial / has integrable multiplicity m 
with respect to a derivation Z), if m is the largest integer for which the following 
is true: there are m — 1 exponential factors exp(gj//^), j = 1, ... ,m — 1, with 
deg (7j < j. deg /, such that each gj is not a multiple of /. 

We say that a Darboux polynomial / of degree N has algebraic multiplicity m with 
respect to a derivation D, if m is the greatest positive integer such that the mth 
power of / divides £jv(D). 

C. Christopher, J. LLibre and J.V. Pereira in |CLP07| show that these multiplic- 
ities are equal when we consider absolutely irreducible Darboux polynomials. This 
is a deep result, but unfortunately nowadays there is no simple characterization of 
exponential factor cxp{g/ f) when / is reducible. If we can characterize exponential 
factors with the ecstatic curve then perhaps we will compute efficiently an integrat- 
ing factor corresponding to a Liouvillian first integral. 



6.2. Inverse integrating factor. An inverse integrating factor is a Darboux poly- 
nomial with cofactor div{A,B). An inverse integrating factor R has the following 
interesting property: The algebraic limit cycles of the polynomial vector field cor- 
responding to D are factors of R, see e.g. |GLV96| . For other results we can read 
e.g. |CGGL03] . We remark easily that: 



R is an inverse integrating factor 




(6.1) ■4=» AdxR + BdyR = div{A,B)R. 

For a given integer N we can compute, if it exists, R with deg R < N. Indeed, 
we just have to solve the linear system (|6.1|) . With this strategy and with classi- 
cal tools of linear algebra we can compute R with 0{N^) arithmetic operations if 
N>d. 



This kind of linear system also appears when we study absolute factorization, 
see |Rup86| |Rup99| IGao03| ICL07| ISchOTj . Indeed, we can compute the absolute 
factorization of a given polynomial R with a solution {A, B) of ()6.1|) . 
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In [CL07j the authors use this kind of hnear system and show that the absolute 
factorization of R can be performed with 0{N*) arithmetic operations. We recah 
that "soft Oh" is used for readabihty in order to hide logarithmic factors in cost 
estimates. Then the question is the following: 

Can we perform the computation of an inverse integrating factor in a deterministic 
way with 0{N*) arithmetic operations instead of 0{N^)1 
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